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The  evolution  of  the  first  and  second  moments  of  the  photon  density 
distribution  for  a two-node  laser  system  was  studied  by  proposing 
appropriate  transition  rate  expressions  and  using  the  Kubo  equations 
developed  for  a Chapman -Kolmogorov  systea.  A numerical  analysis 
technique  was  used  In  solving  the  rate  equations  and  computer  plots 
were  generated  for  various  combinations  of  laser  parameter  values. 
Anomalously  large  fluctuations,  similar  to  those  occurring  In  single- 
mode lasers,  were  observed  during  photon  buildup  of  the  modes,.  Wnen 
only  one  mode  was  oscillating,  a photon  buildup  in  the  non-lasing  mode 
was  seen  to  cause  similar  fluctuations  in  the  lasing  node.  In  all 
cases,  these  fluctuations  were  accompanied  by  a strong  linear 
relationship  between  the  modes. 


PHOTON  STATISTICS  IN  TWO -MODE  LASERS 


I Introduction 


Background 

Since  the  development  of  the  laser  in  I960,  the  statistics  for 
single-mode  lasers  have  been  studied  extensively  using  both  the  semi- 
classical  approach  and  the  fully  quantum-mechanical  theory.  Calcula- 
tions have  been  made  of  the  single-mode  photon  distributions  below,  at, 
and  above  threshold.  Calculations  have  also  been  made  of  the  time- 
behavior  of  the  moments  of  the  photon  distribution  as  the  laser  operates 
away  from  equilibrium.  These  calculations  have  subsequently  been  ver- 
ified by  experiment. 

The  statistics  for  multimode  lasers  are  not  as  well  developed. 

Only  the  two-mode  laser  has  been  studied  in  any  detail}  and  then,  using 
only  the  semiclassical  approach  which  neglects  the  particle  or  quantized 
nature  of  the  light  field  in  favor  of  a classical  electric  field.  The 
fully  quantum-mechanical  treatment  has  just  recently  begun  to  be  analyzed. 
Purpose 

The  purpose  of  this  work  was  to  extend  the  study  of  two-mode  lasers 
by  investigating  the  photon  distributions  in  equilibrium  and  the  time 
behavior  of  the  first  and  second  moments  of  the  distributions  away  from 
equilibrium  for  a two-mode  laser.  Of  particular  interest  was  the  be- 
havior of  the  system  (for  several  pump  levels  and  degrees  of  mode 
interaction)  during  buildup  of  photons  from  a photon  vacuum  and  during 
transition  from  unstable  to  3 table  equilibrium.  Transition  from  various 
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non-equilibrium  points  waa  also  considered. 

Approach 

To  study  the  behavior  of  the  laser  system,  it  was  assumed  that  the 
system  could  be  described  by  the  Chapman-Kolmogorov  master  equation. 
Therefore,  the  resultant  general  evolution  equations  developed  by  Kubo 
were  applicable.  Transition  rates  which  included  terms  to  account  for 
spontaneous  and  stimulated  emission,  as  well  as  mode  interaction  and 
self-saturation,  were  proposed  for  the  laser.  These  were  then  inserted 
into  the  general  evolution  equations,  and  rate  equations  for  the  first 
and  second  moments  of  the  photon  density  distributions  were  obtained. 

Appropriate  values  corresponding  to  various  combinations  of  pump 
levels  and  degrees  of  mode  coupling  were  chosen  for  the  laser  para- 
meters. For  each  combination,  the  pump  was  considered  constant  while  the 
coupling  was  assumed  to  take  place  through  the  mode  population  densities. 
The  equilibrium  points  for  each  set  of  laser  parameters  were  then  ob- 
tained analytically  from  the  rate  equations  and  were  classified  using 
linear  vibration  theory  to  predict  the  behavior  of  the  photon  densities 
around  each  point. 

To  solve  the  moment  rate  equations  dynamically  and  observe  the 
laser  system  through  transition,  a computational  scheme  was  developed 
which  included  a time- incremented  numerical  solution  of  the  equations. 
Finally,  using  a series  of  initial  photon  densities  with  associated 
variances,  this  program  was  run  for  each  combination  of  parameter  values 
and  several  plots  were  obtained. 
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II  Theory 


Background. 

Just  as  there  are  two  basic  approaches  to  quantum  mechanics,  there 
are  two  »in  methods  for  deriving  laser  fluctuations  fully  quantum- 
nechanically.  One  method  (associated  more  with  the  Schroedinger  pic- 
ture) utilizes  a density  operator  to  derive  a "master  equation".  The 
other  method  (associated  more  with  the  Heisenberg  picture)  consists  of 
"adding  Langevin  operators  to  the  Heisenberg  equation  of  the  field 
operator  and  of  the  atomic  operators"  (Ref  9*273)*  The  two  methods 
are  equivalent  in  that  they  solve  the  problem  equally  well.  He  : ever, 
the  master  equation  approach  is  more  appropriate  in  that  the  resulting 
non-linear  equations  can  be  reduced  to  a form  solvable  by  numerical 
methods,  whereas  the  Langevin  method  requires  analytical  solutions 
which  restrict  its  use  to  linear  or  linearized  equations  (Ref  9*273). 
Details  of  both  methods  can  be  found  in  the  work  of  Sargent,  Scully, 
and  Lamb  (Ref  11). 

Single-Mode  Theory.  Using  the  master  equation  approach  for  single- 
mode lasers,  the  probabilities  cf  n photons  existing  are  determined  by 
the  diagonal  elements  of  the  associated  field  density  matrix 

(Ref  Hi 284),  Therefore,  the  equations  of  motion  for  these  elements  det- 
ermine the  flow  of  probability  into  and  out  of  each  n statei 

ijt)  =.  - + 

— V*')  farXO  '/(o  V 

where  oC  is  the  linear  gain  coefficient,  V is  the  loss  coefficient, 
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f3  is  the  saturation  coefficient,  and  the  sign  (positive  or  negative) 
of  each  term  on  the  right-hand  side  indicates  probability  gains  or 
losses,  respectively,  for  the  state  n (Ref  10|2424).  The  photon  dis- 
tribution is  then  determined  by  solving  these  equations  (Eqs  (1)).  It 
is  important  to  note  that  the  diagonal  elements  are  coupled  only 

to  other  diagonal  elements.  Therefore,  to  determine  the  photon  statis- 
tics, the  off-diagonal  elements  of  the  density  matrix  can  be  ignored. 

In  studying  the  behavior  of  any  system,  it  is  often  of  interest  to 
consider  the  behavior  of  the  moments  of  the  distribution  defined  by 

n -o 

where  is  synonymous  with  />Art  above.  When  Eq  (1)  is  substituted 
into  the  time  derivative  of  Eq  (2),  the  result  is  a set  of  equations 
of  motion  for  the  moments  of  the  photon  distribution  which  can  then  be 
solved  numerically  to  observe  the  behavior  of  the  moments  at  as  well  as 
away  from  equilibrium.  Of  particular  interest  are  the  first  and  second 
moments  ( k-1  and  k-2,  respectively)  since  they  can  often  be  deter- 
mined experimentally  for  comparison,  whereas  the  distribution  function 
cannot. 

The  behavior  of  these  moments  has  been  calculated  from  the  initial 
condition  that  no  photons  exist  to  equilibrium  conditions  (Ref  10 j 2426- 
242?).  The  results  for  the  expected  value  of  the  photon  number  (first 
moment)  show  an  initial  slow  rise  dominated  by  spontaneous  emission, 
followed  quickly  by  a rapid,  almost  linear,  increase  as  the  spontaneous 
emission  is  amplified  by  stimulated  emission.  However,  as  equilibrium 
is  approached,  the  saturation  effects  of  the  laser  begin  to  modify  the 
system  gain  and  the  rate  of  increase  of  the  average  photon  number 
approaches  zero.  The  root  mean  sq\*e  deviation  is  determined  by 


combining  the  first  and  second  moments j 


(3) 


S =‘\/<nl>  - 

Thi3  deviation  initially  tends  to  follow  an  exponential  rise  as  for 
thermal  (incoherent)  light  and  exceeds  the  equilibrium  value  as  not  only 
the  photon  numbers,  but  also  the  fluctuations  in  photon  numbers,  are 
amplified.  Then,  as  the  system  moves  into  the  saturation  region,  amp- 
lification of  already  stimulated  photons  becomes  predominant  and  the 
deviations  diminish  to  their  relatively  small  equilibrium  values.  It 
is  important  to  note  that  when  this  problem  was  solved  semiclass ically, 
using  Langevin  forces  to  account  for  noise,  similar  results  were 
obtained  (Ref  9«262-2b3). 

Two-Mode  Theory . Unlike  single-mode  lasers,  the  fully  quantum- 
mechanical  treatment  of  multimode  lasers  has  only  very  recently  begun 
to  be  analyzed  (Ref  8).  Therefore,  for  these  lasers  the  theory  remains 
centered  primarily  around  Lamb's  paper  (Ref  7)  which  gave  the  first 
comprehensive  treatment  of  laser  action.  Following  this  semiclassical 
theory,  the  equations  of  motion  for  a two-mode  laser  system  are 

i,  = ai,(o..-/s.  i.  - ©,aj  (4) 

and 

ix  = (5-) 

where  1^  and  are  the  dimensionless  intensities  of  the  electric  fields 
for  the  two  modes,  a^  and  a2  are  the  gain  coefficients,  andy5*0 
are  the  self -saturation  coefficients,  and  and  are  the  cross- 

coupling coefficients.  When  these  equations  are  integrated  simultan- 
eously for  a range  of  initial  intensities,  the  result  is  a set  of  phase 
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r 


curves  in  the  Ij-Ig  plane  (Refs  7iAl441-Al442;  and  11»125-127). 

• • 

The  possible  equilibrium  points  are  determined  by  setting  1^  - ■ 0 

and  solving  (4)  and  (5)  together.  and/or  I2  equal  to  zero  are 
obvious  solutions,  with  single  frequency  solutions  of 


and 


- (<W/0 
( r, , O * ( a7/.  > o) 


(l) 


(l) 


For  both  modes  to  oscillate  at  equilibrium,  the  parenthesized  expressions 
in  (4)  and  (5)  must  be  zero.  Therefore, 


and 


/s;  x,  +-  olz  = a, 

A T,  +•  I,  = a-t 


(0 


(1) 


which  are  linear  expressions  in  1^  and  I2  and  have  the  solutions 

[ °-|  — (®1.  i]  //S, 


I.  = 


1 -c 


(lo) 


I 


1 - C 


00 


where  the  coupling  constant  is  defined  as 


C = A A , , 

/?,  A (A 

For  this  solution  to  be  a physically  realizeable  point,  both  1^  and  I2 
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Figure  1.  Phase  curves 
showing  the  transient  be- 
havior of  two  mode  oscil- 
lation. Lines  Li  and  L2 
represent  Eqs  (8;  and  (9) 
with  coefficients  aj  - 1, 

■ 1 , ■ 2 , and 

©k-  Qu-  1.  (Ref  11»126) 


Figure  2.  Diagram  similar 
to  Fig.  1,  except  the  gain 
parameter  for  the  second 
mode  is  now  a£  - C.4,  and 
mode  1 inhibits  mode  2 at 
equilibrium.  (Ref  11  * 125) 


must  be  non-negative  values.  When  C <C  1 » this  occurs  when  both 

0-',  - O-,  - (©>.,  40  (13) 

and 

Q-l  = nl  - ) a,  (1 4) 

are  non-negative.  See  Fig.  1 (Ref  11 1 126).  aj  and  a^  are  called  the 
effective  gains  of  their  respective  modes.  If  either  effective  gain 
is  negative,  then  the  associated  mode  will  not  oscillate  and  the  laser 
will  lapse  into  single-mode  operation  as  shown  in  Fig.  2 (Ref  11 1 125) 
with  mode  intensities  as  defined  by  either  (6)  or  (7),  as  appropriate. 

In  this  case,  assuming  both  and  a2  are  positive,  the  oscillating 
mode  is  said  to  inhibit  oscillation  in  the  other  mode. 
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When  C>  1 , the  intersection 
of  the  tvo  lines  (8)  and  (9)  Is  a 
physical  point  only  when  both  aj  and 
a,  are  non-positive j that  is,  vhen 
each  mode  tends  to  inhibit  the  other, 
In  this  case,  the  Intersection  point 
is  an  unstable  equilibrium  and,  as 
with  the  Inhibition  of  one  mode,  the 
system  relaxes  into  single-mode  op- 


eration. Now,  however,  both  single- 
mode solutions  defined  by  (6)  and 
(7)  are  stable  solutions  and  the 


Figure  3.  Phase  curves 
showing  the  transient  be- 
havior of  two-mode  oscil- 
lation with  coefficients 
al  - a?  - 1 , /?,  - /3t-  1 , 
and  ©Ii.  ■ " 2 , (Ref 


point  to  which  the  system  evolves  11 i 127) 

depends  on  the  initial  conditions. 

See  Fig.  3 (Ref  Ilil27). 

When  C - J.  , lines  (8)  and  (9)  are  parallel  and  Eqs  (10)  and  (11) 
are  Indeterminate.  In  this  instance,  any  combination  of  intensities  on 
either  mode  equilibrium  line  (as  defined  by  Eqs  (8)  and  (9))  will  be 
stable.  Along  these  lines,  the  stability  will  be  neutral  with  any  de- 
viations neither  damped  nor  magnified.  Deviations  away  from  these  lines, 
however,  will  decay  in  time,  bringing  the  system  back  to  the  line. 

As  can  be  seen,  the  degree  of  coupling  between  modes,  as  expressed 
by  the  coupling  constant  C,  is  an  important  factor  in  determining  the  be- 
havior of  the  system.  The  critical  point  is  when  C - 1 . Therefore, 
three  regions  of  coupling  strength  are  defined  i 


weak  coupling 
neutral  coupling 


c<  1 

C = 1 


8 


0*) 

(/sL) 


c > 1 


strong  coupling 


r 


(Vc) 


where  the  degree  of  coupling  applied  to  each  region  was  determined  by 
considering  Sq  (12)  and  the  relationship  of  ©,^6^  to  (3,  . Weak 

coupling  implies  that  the  cross-coupling  coefficients,  while  they  exist 
and  affect  the  system,  are  not  the  dominant  terms.  Instead,  the  self- 
saturation coefficients  are  primarily  responsible  for  limiting  the 
growth  of  the  system.  Qualitatively,  this  is  expressed  as 
or  C<  1 . Conversely,  strong  coupling  indicates  that  cross-coupling 
is  the  dominant  effect,  and  is  expressed  as  • Neutral 

coupling,  then,  corresponds  to  neither  effect  dominant,  or  . 

Derivation  of  Moment-Rate  Equations 

General  Evolution  Equations.  The  two-mode  laser  (or  more  correctly, 
the  two  laser  modes)  can  be  considered  as  a system  composed  of  many 
elementary  units  which  interact  with  themselves  and  with  the  environ- 
ment. As  such,  a set  of  macrovariables  is  defined! 


X = ( x,,  O 


(1C) 


where  and  Xo  represent  the  photon  numbers  present  in  modes  1 and  2, 
respectively.  It  is  assumed  that  X is  a stochastic  variable  which 
follows  a Markovian  process  as  it  evolves  in  time.  Therefore,  the 
master  equation  of  the  system  is  the  Chapman-Kolmogorov  equation,  which 
may  be  written  as 

JL.p(x,t)  =-/vJ  (uVfo,t)4i  +jW(K-£,£)Ki-i:.t)'lE  O7) 

where  P(X,t)  is  the  probability  distribution  of  X (i.e.,  the  probability 
that  there  exist  X^  photons  in  mode  1 and  X2  photons  in  mode  2)  at  time 


r 


t,  and  V(X,r)  Is  the  transition  rate  for  a jump  froa  X by  an  arount  r. 
The  evolution  of  variables  governed  by  this  equation  has  been  studied 
extensively,  beginning  with  the  work  of  van  Kaapen  (Ref  16)  and  contin- 
ued by  Kubo  (Refs  iJ-j  5j  and  6)  and  Haken  (Ref  3),  Therefore,  the 
resulting  evolution  equations  for  the  systea  are  just  listed  herei 


i(t)  = £,(4.,^ 


(nr) 


and 


i(t)  +-  rv3cl(y,t') 


+ 


(n) 


where /)_  is  the  slae  of  the  systea,  x(t)  is  the  expected  value  of  the 
normalized  intensive  aacrovariable  (representing  the  photon  densities) 
of  the  systea,  *(tr)  - X / CL  , *nd  (t)  is  the  covariance  tensor 
for  x.  The  variables  and  c2  , as  defined  by 


C.,  ( i w(^)rN)£ 


(20) 


r r 


UO 


are  the  first  and  second  aoaents  of  w(Irr),  the  transition  rate  for  the 
intensive  variable  xi 


oO 


- ^1,0  - W ( X > £)/sL 


(22) 


1 
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, A derivation  of  Eqs  (16)  and  (19)*  starting  froa  (17)  and  parallelling 

the  single  variable  derivation  of  Kubo  (Ref  5),  is  included  in  Appendix  A. 

A general  analysis  of  Eqs  (18)  and  (19)  predicts  sone  interesting 
behavior  on  the  part  of  the  laser  system.  First,  the  expected  values 
of  the  photon  number  densities  follow  deterministic  paths.  For  a given 
set  of  laser  operating  conditions,  the  path  to  equilibrium  is  determined 
solely  by  the  initial  photon  densities.  For  a non-linear  relaxation, 
the  fluctuations  generally  exhibit  unusual  behavior  as  the  system  moves 
toward  equilibrium.  This  effect  becomes  very  pronounced  when  the  system 
moves  froa  the  vicinity  of  an  unstable  equilibrium  toward  a stable  eq- 
uilibrium. Around  each  point*  both  the  photon  number  densities  and 
their  fluctuations  grow  or  decay  exponentially,  depending  on  tne  stability 
of  the  point.  In  between,  the  fluctuations  grow  anomalously  large,  peak- 
ing about  half-way  to  the  final  equilibrium.  (Refs  6i?8-79>  and  4»6 -8) 

This  enhanced  fluctuation  is  predicted  for  and  has  been  observed  in 
single-mode  lasers  (Ref  2).  It  has  also  been  predicted  for  two-mode 
lasers  by  Arecchi  and  Ricca  in  a paper  which  used  a generalized  version 
of  the  Vol terra -Wiener  functional  technique  to  study  multimode  laser 
action  (Ref  l). 

Transition  Rates . To  make  the  evolution  equations  (18)  and  (19) 
specific  to  the  two-mode  laser  problem  requires  the  proper  specification 
of  the  transition  rates  w(^,r).  It  is  important  to  note  that  the  size 
of  the  transition,  r,  is  in  terms  of  the  extensive  variable  X or  the 
total  number  of  photons,  whereas  the  dependence  on  the  state  of  the 
system  is  in  terms  of  the  Intensive  variable  x or  the  photon  density. 

The  basic  problem  is  depicted  in  Fig.  4 for  a system  in  the  state 
x*(n^,n2).  a^  and  b^  are  the  upward  and  downward  transition  rates, 

kL 
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Figure  4.  Transition  rates  associated  with  a laser  system  In 
the  Internal  state  x ■ (nj,^).  Upward  transitions  (f)  correspond 
to  photon  ealsslons  while  downward  transitions  (|)  signify 
photon  losses, 

respectively,  for  x1  to  sake  a jump  r1  , and  c and  d are  cor- 
x i &2  n2 

responding  terns  for  x^.  It  Is  assumed  that  photons  are  added  to  or 
deleted  from  the  system  one  at  a time.  Therefore,  the  only  non-zero 
values  for  r^  and  rg  with  associated  non-zero  rates  are  1 and  -1. 
Furthermore,  even  for  these  values,  the  associated  rates  for  one  mode 
will  be  non -zero  only  If  r^  ■ 0 for  the  other  mode.  Therefore,  there 
are  only  four  events  with  non-zero  associated  transition  rates i 
Increases  by  one,  r"(l,0)  ; X^  decreases  by  one,  r ■ (-1,0)  j Xg  In- 
creases by  one,  r-(0,l)  j and  X£  decreases  by  one,  r*(0,-l)  , Based 
on  this  and  on  Fig.  4,  the  transition  rates  for  can  be 

written  as 


These  rates  are  also  valid  for  w(jr,r)  since  jr(t)  is  just  the  expected 
value  of  x. 


* 

> 


* 


How,  then,  should  a , b , c , and  d be  characterized?  b 
"l  "l  “2  "2  °1 

and  d essentially  are  the  loss  terse  of  the  system,  including  internal 

cavity  losses  as  well  as  the  amount  coupled  out  of  the  system.  There- 
fore, they  can  be  expressed  as 


and 


k,  - K n, 

\ - k<v 


Us) 


where  hj  and  h2  are  determined  by  the  cavity  lifetime  for  a photon  in 

■odes  1 and  2,  respectively  (Ref  12i4l3).  a and  c , representing  the 

nl  n2 

gain  of  each  mode  of  the  system,  are  somewhat  more  complicated.  Each 
gain  expression  must  Include  terms  to  account  for  both  stimulated  and 
spontaneous  emission  as  well  as  for  self -saturation  and  cross-coupling. 
The  self -saturation  occurs  because  there  is  only  a finite  pump  which 
excites  a finite  number  of  atoms  or  molecules  which  may  then  be  stimu- 
lated to  contribute  to  the  system  of  photons.  Cross -coupling  occurs 
when  both  modes  have  to  compete  for  some  of  the  sane  atoms  or  molecules. 

The  expression  chosen  which  meets  these  requirements  for  a is 

nl 


a 


3m 

I +■  <*,,0,  i- 


where  g^  and  g21  are  factors  accounting  for  spontaneous  and  stimulated 
emission,  respectively,  and  c*M  and  account  for  self -saturation 
and  cross -coupling,  respectively.  A similar  expression  for  c is 


¥ 


^ + Su.n,. 

I +■ 


with  g12,  f^22*  *»-«  &nd  tevln«  th®  SABe  Beanln«  mode  2 as  g^, 

g21,  c^(1_  and  o^jj.  do  for  mode  l.  It  Is  important  to  note  that  with 
these  expressions  the  transition  rates  will  always  be  non-negative.  The 
validity  of  these  choices  will  be  demonstrated  in  a later  section. 

Laser  Rate  Equations,  When  the  transition  rates  (Eqs  (24) -(27)) 
are  inserted  into  the  evolution  equations  (18)  and  (19)*  the  resultant 
equations  are  the  desired  laser  rate  equations  for  the  first  and  second 
moments  of  the  photon  density  distributions.  For  ease  of  manipulation, 
the  evolution  equations  will  be  handled  in  their  component  fora.  There- 
fore, Eq  (18)  becomes 


ji  W = c-.i 


and  Eq  (19)  is  written  as 


— C 

~SL  Sk 


+ 


where  now  c^(^)  and  c2(^)  are  expressed  as 


I 

and 

" for 


Uo) 


Dealing  first  with  Eq  (30)*  Eqs  (24)-(2?)  are  appropriately  inserted 
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into  Eq  (30)  as  determined  "by  Eqs  (23).  The  resulting  components  of 
Cj^jr)  are 


C 


U 


(-Ob* 


r h—+.  n-i - \^x<\^  (32.) 

I +■  0^ 

and 

e-xAl)  - fkcc^C^O 

~ C ^ C ~ ^ ^ 

- ii2-.±  (33) 


Similarly,  the  diagonal  elements  of  c ^ (i)  are  found  from  Eq  (31)  to  be 


c = 3"  ^ A'  + b.rt, 

" l + WIV(\U 

and 

C_  _ ^ fl*- 4-  'n2_'n1__ 

Ztl-  \ -V  t-  o^nO,. 


Since  mutually  exclusive  events  were  assumed,  r^  or  r^  must  equal  zero 
whenever  j f k . Therefore,  the  off-diagonal  elements  of  c ^ are 
both  zero. 
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To  complete  the  development  of  the  rate  equations,  the  partial 
derivatives  of  (32)  and  (33)  with  respect  to  y^  and  y2  must  be  deter- 
mined. In  doing  this,  as  throughout  this  development,  it  should  be 
pointed  out  that  n^  and  it,  are  defined  as  synonymous  with  yj  and  y2, 
respectively.  When  this  is  done  and  Eqs  (32)-(35)  are  Inserted  into 
Eqs  (28)  and  (29)  as  appropriate,  the  rate  equations  for  the  first  and 
second  moments  of  the  photon  density  distributions  become 


A 


\ 


V 4-  Si,  A, 
\ 4- 


+ 5i-i.ni. 
I +•  Av 


r — 

/ \~ 

°Al  +■  !)il  A L ^ 

( l ' 

(i+<Al+«<aniy 

. — L/  Su  +V.A)  ( 35r) 

JT  l lvo/hn,  / 


a 

LV  " +■  Ai  ^ 

- <2^ 

, 1 w4l.r\Ly 

In 

1 

m 

+ 3*i. 

1+  A,d- 
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and 


^ = -*;» 


/ \” 

p* 

/ N — 

+ Vi  ) 

rr 

°^>2_  \ ^ ^ 2.1  ) 

0 + 

( U*„n, 

+ <rll 


( i +■  1 *• 


There  Is  actually  a sixth  rate  equation,  for  <3^  , hut  It  Is  Identical 
In  fora  to  (40).  Also,  the  equality  of  and  0^  due  to  the  sym- 

metry of  the  covariance  matrix  is  used  to  combine  terms  in  Eqs  (38), 
(39),  and  (40).  Therefore,  only  the  five  equations  (36)-(40)  need  be 
solved  simultaneously  to  determine  the  behavior  of  the  laser  system. 


Equilibrium  Analys is  of  Moment  Equations 

As  stated  previously,  the  purpose  of  this  study  was  to  determine 
the  behavior  of  a two-mode  laser  system  as  it  transitioned  from  non- 
equilibria to  equilibrium.  This  was  the  reason  for  developing  the 
laser  rate  equations  and  solving  them  dynamically.  However,  before 
applying  numerical  methods  to  these  equations  to  follow  the  system 
through  its  transition,  a basic  understanding  of  the  system  can  be  gained 
by  determining  analytically  the  positions  of  all  equilibria,  both  stable 
and  unstable,  and  the  behavior  of  the  system  around  each  of  these  points. 
Of  particular  Interest  is  the  first-moment  equilibria,  since  these 
values  for  n^  and  Tir,  are  necessary  to  do  any  equilibrium  analysis  on  the 
variances.  Therefore,  the  photon  rate  equations  (36)  and  (3?)  are  con- 
sidered first,  followed  by  a brief  discussion  of  the  variance  rate 
equations  (38),  (39),  and  (40)  in  equilibrium. 

First-Moment  Equations.  Recalling  the  first-moment  or  photon 
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density  rate  equations  (36)  and  (37).  equilibria  will  occur  when  both 
and  are  zero.  These  are  also  referred  to  as  singular  points  of  the 
system.  The  point  n^  ■ ■ 0 can  only  be  a solution  if  g^  ■ g^2  “ 0 j 

but  this  corresponds  to  zero  noise  or  spontaneous  emission  in  the  system, 
which  is  an  unrealizeable  situation.  In  fact,  neither  n^  nor  ^ equal 
tp  zero  can  be  a solution  for  finite  values  of  n.  Therefore,  there  is 
some  combination  of  non-zero  (n^,^)  for  which  equilibrium  will  exist. 

As  a first  approximation,  assume  n^  and  1 , g^  - g21  , and 

g12  " g22  * Then  51(13  (36)  and  (37),  set  to  equilibrium,  become 


A,  — o = Y\x 


(4  0 


Aj.  — 0 - “h,.  Oz  (42.) 

Since  only  non-zero  values  of  n^  and  are  of  interest,  the  bracketed 
quantities  in  (41)  and  (42)  are  set  to  zero  and  examined.  After  re- 


arranging, the  resultant  equations  are 


4-  - I 


(45) 


= * 


6*4) 


These  are  Just  two  straight  lines  (call  them  1 ^ and  12,  respectively) 
in  (n^,!^)  space  with  r^-  and  n^-  intercepts  of 
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% 

* 


for  1^,  and 


for  lg. 

Within  the  limits  of  the  original  assumption  (n1,n2»l),  the 
point  of  intersection  of  these  two  lines  is  the  location  of  one  equil- 
ibrium or  singular  point.  To  find  this  point  of  intersection,  equations 
(43)  and  (44)  are  solved  simultaneously.  This  results  in 


Now,  if  the  following  substitutions  are  made, 


1 

H 

f> 

~ 

(41  *) 

(41b) 
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K ^ 


equations  (47)  and  (4ti)  become 


n,  = 


k-(e.,/iOO/A 

I - C. 


and 


/l. 


I - c 


where 


c - 


°<|».  Q<».| 


u. 


(3,  fi. 


(41c.) 


(iTo) 


< 5"/  ) 


(^<0 


is  a measure  of  the  degree  of  coupling,  as  previously  discussed.  Com- 
paring this  to  the  semlclassical  development  of  Lamb  presented  earlier, 
it  is  seen  that  Sqs  (50)  and  (51)  ere  essentially  identical  to  Eqs  (10) 
and  (11)  with  the  normalized  intensities  1^  and  I2  replaced  by  their 
quantum  equivalent,  the  photon  number  densities.  Furthermore,  when 
Eqs  (41)  and  (42)  are  expressed  as  single  fractions  (with  non-zero 
denominators),  they  can  be  reduced  to 


and 


» 0 * ", 


A t - 0 - 


(%”0  "k^n,  nt 
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which  are  equivalent  to  Lamb's  equations  of  motion  for  1^  and  , Eqs 

(4)  and  (5)»  at  equilibrium.  Therefore,  in  the  limits  of  equilibrium 

and  of  large  n (or  zero  noise),  the  chosen  transition  rates  an  * t 

c , and  d are  valid,  and  the  semiclasslcal  equilibrium  analysis  done 
"2  °2 

by  Lamb  can  be  applied. 

When  the  spontaneous  emission  terms  and  g^  are  not  assumed  away, 
the  singularities  for  the  system  are  determined  by  setting  Eqs  (36)  and 
(37)  to  zero  and  solving  simultaneously.  Solving  Eq  (36)  for  n 2 and 
substituting  the  obtained  expression  into  (37),  the  result  is  the  fourth 
order  equation 


+■  ~v»,y 


^ - 0 (sy) 


This  is  then  solved  using  basic  algebraic  techniques  (Ref  14il07)  and, 
depending  on  the  coefficients,  generally  results  in  four  singular 
points  for  the  system. 

Characterization  of  Singular  Points.  The  behavior  of  the  system 
around  a given  singularity  is  determined  by  considering  a small  dis- 
turbance around  that  point.  Based  on  this  behavior,  the  singularity 
can  then  be  classified  into  one  of  six  basic  types,  as  shown  in  Fig.  5. 
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equations  can  be  combined  in  the  fora 


dn, 

Q-  -+■  bn, 

(sx) 

<K 

c.  n'  +-  d nt' 

where 

a = -h|0^ 
b = ( ~ ^ > 
q-  = (wk 

K / ( \ +■  <*..  ^1°  + 0°  ^ 

- zln,o<u  n,  - ln,^i7.^l^//j  + o<'unic>+- 

- ZV^^n*-  n,*)/ V 1 +«^*i of  Wnnt°) 

(nO 

(s'?  o') 

and 


l + o,°  + cx'tt  nL*N) 


The  type  of  singularity  can  then  be  determined  by  a series  of  expressions 
involving  these  coefficients.  These  expressions,  and  their  relation- 
ships to  the  types  of  singular  points  as  shown  in  Fig,  5,  are  as 
follows  (Ref  15*4^4-)  t 


I. 

(k 

-0 

Z+4ad  > 0 

f ad  - b c < (9 

|^ad  - be  > 0 
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(k 
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Mode.  ( Ft^*  £c) 

Also, 

a node  or  spiral  is  unstable  or 

• stable  depending  on 

Nlode 

Saddle.  ( Fij.  5d^ 


is  positive  or  negative,  respectively. 
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Second -Moment  Equations.  Knowing  the  equilibrium  values  for  the 
photon  densities,  the  variance  xate  equations  (38)#  (39)#  and  (40)  can 

e 

be  set  to  equilibrium  and  solved  simultaneously.  The  zero  of  in 

Sq  (38)  represents  a plane  parallel  to  the  0^  axis  in  ( (Tl,  , (T^ , ) 

space.  Similarly,  Eq  (39)  represents  a plane  parallel  to  the  CTJJ  axis, 
whereas  Eq  (40)  is  the  equation  for  a straight  line.  Therefore,  for  a 
given  singular  point  (n^.n^),  there  will  be  at  most  a single  equilibrium 
point  for  the  variances.  Eqs  (38)  and  (39)  are  solved  for  CT^  and  CT^_ , 
respectively,  to  get 


_ -(iBau-l-c) 

" zfl 

C6  0) 

and 

_ -(lE  <r„.  +C) 
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and. 


~ -H-  .(  1 0,  *+-  V>tnL  (^^4') 


Sqs  (oO)  and  (6l)  are  then  Inserted  into  Eq  (40)  and  the  resulting 
expression  for  is 


CEb  ->-  flftF 


(a  -e 


(63) 


This  value  is  then  used  in  Eqs  (60)  and  (61)  to  obtain  corresponding 
values  of  0"[  and  (T^  , respectively. 
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Ill  Approach 


Computational  Scheme 

In  studying  the  behavior  of  the  laser  system  using  the  laser  rate 
equations  (36) -(40),  several  steps  were  accomplished.  The  first  of 
these  was  the  development  of  a computer  program  which  Included  the  selec- 
tion of  a numerical  analysis  scheme  capable  of  solving  the  system  of 
rate  equations  (36) -(40)  in  arbitrary  time  increments.  The  routine 
chosen  for  this  most  important  procedure  is  a FORTRAN  IV  subroutine 
called  BLCKDQ  (Ref  13).  It  uses  a predictor-corrector  formula  of  8th 
order  with  user  specified  tolerances,  which  were  set  to  10  because 
of  initial  rapid  changes  in  the  moments.  Time-incremented  solutions  are 
possible  because  BLCKDQ  solves  the  system  of  equations  at  a user  spec- 
ified time  (XOUT),  given  appropriate  initial  values  at  some  prior  user 
specified  time  (XIN).  Therefore,  the  system  is  incremented  to  equil- 
ibrium by  calling  BLCKDQ  within  a loop  where  the  "initial"  values  are 
set  equal  to  the  solutions  from  the  previous  time  step,  XIN  is  set  equal 
to  the  previous  XOUT,  and  XOUT  is  incremented  by  some  predetermined  amount. 
Built  around  this  loop,  the  remainder  of  the  program  was  then  developed 
to  solve  the  rate  equations  to  equilibrium  for  a series  of  initial  value 
points  and  to  store  for  plotting  all  sets  of  solutions  belonging  to  a 
given  set  of  laser  parameters.  The  laser  rate  equations  are  contained 
in  a user  supplied  subroutine  FCN  which  is  called  by  BLCKDQ. 


Selection  of  Values  for  Laser  Parameters 

The  second  step  in  making  the  rate  equations  solvable  was  to  deter- 
mine appropriate  values  for  the  various  laser  parameters  (g  .,  , h. ) 

for  the  desired  combinatins  of  pump  levels  and  degrees  of  coupling. 
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To  do  this,  several  assumptions  were  made.  First,  it  was  assumed  that 
the  cavity  losses,  the  self -saturation,  and  the  cross -coupling  were  the 
same  for  both  modes.  That  is 


u = L 

( ) 

(64  b) 

*,a  = ^a, 

(64e) 

Comparing  Siegnan’s  photon  rate  equation  (Ref  12,418) 

to  Sqs  (36)  and 

(37),  and  noting  that  Siegman  has  no  saturation  or  coupling  terms,  the 

following  equivalences  can  be  mde, 

3 a,  “ K,(nIu  ' Mj^ 

1>  a) 

^ a.  ( isl  ^ - rOjL  ^ 

> 

(lS  0 

(( e ^ c ) 

= K t k1u 

(6^=0 

where  and  are  population  densities  of  the  upper  and  lower  lasing 
levels  and  and  are  coupling  coefficients  between  the  lasing  modes 
and  the  lasing  medium.  It  is  assumed  that  the  lower  lasing  level  de- 
populates essentially  Instantaneously  ( s;  0 ) and  the  upper  lasing 
level  is  maintained  constant.  Therefore,  - g21  - constant  and 

g12  “ g22  " cons'tan't* 

To  select  a value  for  h,  the  following  relationship  is  noted 
(Ref  12,418) j 

V\  - i/-re  l'u'> 
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where  is  the  photon  cavity  lifetime,  A representative  value  for  'XCm 
S 

is  10  sec.  Therefore 


(41) 


is  used  for  all  combinations  of  g^  and  , The  saturation  and 

coupling  coefficients  are  somewhat  more  elusive.  To  get  an  idea  of  at 
least  an  order  of  magnitude,  the  expressions  (49)  were  recalled  and 
compared  with  values  for  ^ and  O ^ used  by  Lamb  (Ref  lit  121). 
Therefore,  for  each  type  of  coupling  (as  determined  by  (6)  and  (52)), 
the  following  values  were  chosen: 


weak  coupling  ( C.  <C  1 ) » <X  „ - - 2.  x \0 

= lO-,° 


- NO 


neutral  coupling  (c*l)»  -10 

i 


- to 


strong  coupling  ( c_  > 1 ) t ~ °4»-  " 


.10 


o<(t  - <,  - 2.k*0 


-(O 


(41) 

(70) 


In  selecting  appropriate  pump  parameters  g 
levels  were  defined: 


ij  ’ 


five  different  pump 


well  above  threshold  : 
slightly  above  threshold  : 

threshold  : 
slightly  below  threshold  : 
well  below  threshold  : 


3*v>  2_Vn  x 

K<  ^t(< 

^ = v*, 

o*^,< 5-.  < V), 

3*.  <O.S^ 


^ ».  ^ «•«.  K.  f*  I 
^ <o.sV 


For  each  pump  level,  both  g21  - g22  and  g21  > g22  were  considered. 
For  the  case  of  g2^  > g22  , the  values  were  chosen  so  that  mode  2 was 
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totally  Inhibited  by  mode  1.  This  corresponds  to  n,  being  negative  at 
the  equilibrium  defined  by  (48).  The  values  chosen  for  the  pump  para- 
meters are  found  in  Table  I,  as  veil  as  the  other  laser  parameter  values 
previously  mentioned. 

The  only  other  variable  in  the  laser  rate  equations  which  must  be 
specified  is  the  size  of  the  system, -fi_  . The  only  restriction  on  _0_ 
is  that  it  be  some  finite,  positive  value.  Therefore,  for  convenience, 
a size  of  _TL=1  was  chosen. 

Selection  of  Initial  Conditions 

For  each  set  of  laser  parameters  the  final  step  before  running  the 
program  and  plotting  the  results  was  to  determine  a series  of  initial 
value  points  in  ( n^ , ^ ) and  their  associated  variances.  The  general 
goal  was  to  select  between  20  and  25  starting  points  per  set  of  para- 
meters spread  around  the  n^  vs,  n^  plot  area.  To  do  this  required  the 
determination  of  a maximum  photon  number  NMAX  which  then  became  the 
dimension  (in  (n^iig)  space)  for  the  square  plotting  area.  For  cases 
with  pump  levels  above  threshold,  NMAX  was  taken  to  be  the  maximum  inter- 
cept value  of  1^  and  1^  as  defined  by  (45)  and  (46).  For  threshold  and 
below  threshold  cases,  NMAX  was  determined  by  identifying  the  stable 
equilibrium  for  the  particular  case  and  multiplying  the  larger  n-  coor- 
dinate by  an  arbitrary  factor  between  approximately  three  and  ten.  This 
ensured  that  the  stable  equilibrium  was  within  the  plotting  area,  but 
was  not  obscured  by  the  axes.  The  values  chosen  for  NMAX,  as  shown  in 
Table  I,  were  constant  for  each  pump  level,  regardless  of  the  degree  of 
coupling.  Initial  points  were  then  selected  evenly  around  the  edge  of 
the  plotting  area  in  addition  to  several  points  around  the  origin  and 
points  close  to  any  pre -determined  unstable  equilibria. 
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Parameter  Value 


matched  modes 

identification  by  which  each  set  of  conditions  will  referred  in  the  plots 


Having  selected  the  initial  (nj,n2)  , the  question  is  how  to  deter- 
mine, or  choose,  appropriate  initial  variances  for  each  point.  To  do 
this,  another  question  was  asked.  How  would  each  initial  point  be 
attained  physically  so  that  the  system  could  then  make  a transition  to 
the  desired  stable  equilibrium?  One  possibility  is  to  first  change 
some  of  the  laser  parameters  such  that  the  initial  point  becomes  a 
stable  equilibrium  point.  The  system  is  then  returned  "Instantaneously" 
to  the  original  configuration  and  the  behavior  of  the  system  in  transi- 
tion is  observed.  This  scheme  was  used,  varying  the  cavity  loss  terms 
h^  and  h2  , as  in  a Q-switched  laser.  The  necessary  "new"  loss  terms 
were  determined  by  solving  Eqs  (36)  and  (37)  at  equilibrium  for  h^  and 
h2  , respectively,  given  and  . Then,  using  these  "new"  values  and 
the  initial  n^  and  n2  , the  associated  variances  were  obtained  by  solv- 
ing Eqs  (63),  (60),  and  (61),  In  several  cases  (above  threshold, 
strongly-coupled),  some  of  the  initial  points  could  not  be  changed  into 
stable  equilibria  by  merely  changing  the  cavity  losses.  In  these  in- 
stances, the  initial  values  for  the  variances  were  arbitrarily  set 
as  follows  1 


(7l0 

= K 

(nO 

(7lc) 

In  running  the  computer  program  to  determine  the  behavior  of  the 
laser  system,  it  was  decided  that  the  standard  deviations  Sj^  and  s2  and 
the  correlation  coefficient  defined  by 
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and 


<r„. 

V ^ ^ ' 

would  be  of  more  Interest  than  the  actual  variances.  Therefore,  these 
values  were  calculated  and  printed  out  along  with  the  photon  densities, 
the  variances,  and  XOUT  each  time  the  program  cycled  through  the  3LCKDQ 
loop.  In  addition,  besides  the  n^  vs.  rig  plots  similar  to  Lamb's 
results,  three-dimensional  plots  Involving  s^,  s^,  and were 
produced . 
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IV  Results  and  Discussion 

Equilibrium  Analysis 

Photon  Densities.  The  photon  density  equilibrium  points  for  each 
combination  of  laser  parameter  values  (as  defined  In  Table  I)  are 
listed  in  Table  II.  Included  Is  the  classification  of  each  point. 

As  expected,  there  are  generally  four  such  points  for  each  set.  The  only 
exception  Is  when  the  modes  are  neutrally  coupled  (i.e.,  when  the 
cross -coupling  terms  equal  the  self -saturation  terms).  In  this  case, 
there  are  either  two  or  three  singular  points,  depending  on  whether 
the  pump  parameters  are  equal  or  unequal,  respectively.  For  all  but 
the  threshold  cases,  these  points  are  very  close  to  the  equilibria  of 
the  noiseless  approximation  (as  defined  by  Eqs  (41)  and  (42)).  The 
difference,  though  small,  is  important  because  it  indicates  that  a 
photon  vacuum  is  never  an  equilibrium,  stable  or  unstable,  for  either 
mode.  Therefore,  even  when  there  is  no  net  gain  for  the  system,  there 
is  some  positive  average-photon -density  to  which  the  system  will  relax. 

For  each  set  of  parameters,  there  is  at  least  one  stable  and  one 
unstable  node.  The  type  of  the  other  singularities  depends  on  both 
the  degree  of  coupling  and  the  level  of  pump.  For  all  weak  and  neu- 
tral coupled  cases,  these  other  equilibria  are  all  saddle  points,  each 
of  which  has  one  negative  coordinate  thus  ensuring  that  the  system  can 
move  toward  stability  from  any  physically  realizeable  point.  When  the 
laser  modes  are  strongly  coupled,  the  type  of  these  other  singularities 
is  determined  by  the  level  of  pumping.  At  threshold,  these  two  equil- 
ibria are  both  complex-valued.  Below  threshold,  both  singularities 
are  unstable,  one  being  a saddle  point  and  the  other  an  unstable  node 
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Table  II.  Classification  of  Photon  Density  Singular  points 


Data  Set 

ID  n, 

n2 

Type 

All 

1 .333x10 10 

I.333xl010 

stable  node 

2.000x1010 

-2.500 

saddle 

-2.500 

2.000x1010 

saddle 

-1.250 

-1.250 

unstable  node 

A12 

2.000X1010 

5.000 

stable  node 

-1.538 

7. 500x1 09 

saddle 

2.l67x1010 

-3.333x109 

saddle 

-1.250 

-1.667 

unstable  node 

A13 

4.000x1010 

1 . 250 

stable  node 

1.333x1010 

1.333xl010 

saddle 

1.250 

4.000x1010 

stable  node 

-1.250 

-1.250 

unstable  node 

A1 4 

4.000x1010 

0.3846 

stable  node 

-5.000 

1 .500x1010 

saddle 

-3.333x109 

2. 167x1010 

stable  spiral 

-1.250 

-1.667 

unstable  node 

A15 

4.000x1010 

1 .000 

stable  node 

-2.003 

1 .500x1010 

saddle 

-1.250 

-1.667 

unstable  node 

A21 

4.000x103 

-27.000 

saddle 

-27.000 

4.000x10S 

saddle 

2.667x10s 

2.667x10S 

stable  node 

-13.500 

-13.500 

unstable  node 

A22 

4.000xl0S 

51.000 

stable  node 

-15.429 

1.000x10s 

saddle 

4.667x10s 

-1.333x10s 

saddle 

-13.500 

-51.000 

unstable  node 
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Table  II.  Classification  of  Photon  Density  Singular  Points  (cont) 


\ 


Data  Set 

ID  n. 

n-, 

Tyne 

A23 

8.000x10® 

13.500 

2.667x10® 

-13.500 

13.500 

8.000xlOS 

2. 667x10S 
-13.500 

stable  node 

stable  node 

saddle 

unstable  node 

A24 

8.000x10S 

-27.000 

-1.333x10S 

-13.500 

7.286 
2.000x10S 
4.667x10S 
-51 .000 

stable  node 

saddle 

stable  spiral 

unstable  node 

A25 

8.000x10S 

-18.000 

-13.5000 

17.000 
2.000X108 
-51 .000 

stable  node 

saddle 

unstable  node 

A31 

1 .OOOxlO5 
-I.OOOxlO5 

5. 7735x1 O4 
-5. 7735x1 04 

-1 .OOOxlO5 

1 .OOOxlO5 
5.7735x104 

5. 7735x1 04 

saddle 

saddle 

stable  node 

unstable  node 

A33 

5.7735X104 

-5.7735X104 

I.OOOxlO5! 

-1 .000x1 05i 

5. 7735x1 O4 
-5.7735x104 
-1 .OOOxlO5! 

1 .OOOxlO5! 

stable  node 

unstable  node 

complex 

complex 

A3  5 7.071 IxlO4  7.071 IxlO4  stable  node 

-7.071 IxlO4  7.071 IxlO4  unstable  node 


A41 

49.000 
-6. 667x1 O7 
-1.000x10s 

98.000 

49.000 

-6.6b7x107 

98.000 

-I.OOOxlO8 

stable  node 

unstable  node 

saddle 

saddle 

A42 

49.000 

11.500 

stable  node 

-49.000 

-4.000x10® 

unstable  node 

-I.OOOxlO8 

13.143 

saddle 

1.333x10S 

-4.667x10s 

saddle 

J 
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Table  II.  Classification  of  Photon  Density  Singular  Points  (cont) 


Data  Set  ID 


49.000 

-6.667x107 

-2.000x10® 

-49.000 


49.000 

-2.000x10® 

-7.000 

-4.667x10® 


49.000 

-2.000x10® 

-16.333 


1 .000 

-1.667x109 

-2.500X109 

2.000 


1.500 
-11.994 
-2. 000x1 09 
3.333x10® 


1 .000 

-1 .667x109 
-5.000X109 
-1 . 000x1 09 


1.500 

-4. 000x1 O9 
-4. 667x1 O9 
-0.4286 


1.500 

-4.000x10^ 

-1.200 


49.000 

-6.667x107 

-49.000 

-2.000x10® 


11.500 

23.000 

-3.000x10® 

1.333x10® 


1 1 . 5000 
15.333 
-8.000x10® 


1 .000 

-1.667x109 

2.000 

-2. 500x1 O9 


0.1111 
-4. 500x1 O9 
0.1429 
-4.667x109 


1 .000 

-1.667x109 

-1.000 

-5.000x109 


0.1111 
1 .000 
3.333x10® 
-9. OOOxl O9 


0.1111 

0.2000 

-9.000x109 


stable  node 
saddle 

unstable  node 
unstable  node 


stable  node 
saddle 

unstable  node 
unstable  spiral 


stable  node 
saddle 

unstable  node 


stable  node 
unstable  node 
saddle 
saddle 


stable  node 
unstable  node 
saddle 
saddle 


stable  node 
saddle 

unstable  node 
unstable  node 


stable  node 
saddle 

unstable  spiral 
unstable  node 


stable  node 
saddle 

unstable  node 
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or  an  unstable  spiral.  Above  threshold,  when  the  two  modes  are  pumped 
equally,  there  Is  a second  stable  node  which  is  physically  realizeable. 
The  fourth  singularity  is  a saddle  point  lying  essentially  at  the  inter- 
section of  lines  1^  and  1^  (Eqs  (43)  and  (44)).  Given  the  initial  con- 
ditions, it  determines  to  which  stable  node  the  system  will  relax.  When 
only  one  mode  is  inhibited,  the  second  stable  equilibrium  point  moves 
to  a non-physical  point  and  becomes  a stable  spiral.  The  saddle  moves 
to  a point  just  outside  the  physical  region  and  maintains  its  relative 
position  between  the  two  stable  points  to  prevent  physical  points  from 
relaxing  towards  the  non-physical  stable  point. 

Second  Moments . For  each  photon  density  singularity,  there  is  one 
set  of  equilibrium  values  for  the  second  moments  (as  determined  by  Eqs 
(60),  (61),  and  (63)).  For  the  unstable  photon-density-equilibria,  the 
associated  second-moment  equilibria  are  non-physical;  at  least  one  of  the 
variances,  or /°22  ' is  negative.  However,  since  the  points  are  un- 
stable in  the  photon  densities,  the  uncertainties  and  thus  the  second 
moments  will  grow  or  diverge  in  time,  and  the  system  will  move  away  from 
instead  of  towards  these  non-physical  variances. 

For  the  stable  photon  singularities,  the  associated  second-moment 
equilibria  are  real,  physically  meaningful  points.  As  such,  they  are 
listed  in  Table  III  along  with  the  associated  relative  standard  devia- 
tions and  correlation  coefficients.  At  and  below  threshold  the  system 
experiences  large  uncertainties  as  the  standard  deviations  are  of  the 
same  order  as  the  photon  number  densities.  These  large  uncertainties 
are  also  experienced  by  modes  above  threshold  which  are  inhibited  by 
the  other  mode  from  lasing.  Once  a mode  begins  lasing,  the  relative 
standard  deviation  at  equilibrium  drops  sharply,  decreasing  as  the 
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photon  number  density  increases.  It  is  interesting  to  note  that  for  a 
weakly-coupled  system  above  threshold,  matched-mode  oscillation  ( n^  * ) 
shows  for  each  of  the  modes  a greater  variance  than  for  the  single  mode 
that  oscillates  and  inhibits  the  other.  This  is  even  though,  for  the 
same  pump  level,  the  photon  density  is  greater  for  the  single  mode  than 
for  either  of  the  matched  modes.  However,  this  general  phenomenon  should 
be  expected  because  of  the  relative  influences  of  the  second  mode;  in 
one  case,  it  oscillates  at  equal  Intensity  while  in  the  other  Instance 
it  exists  barely  above  photon  vacuum. 

The  correlation  coefficient  indicates  a significant  linear 

relationship  between  the  photon  densities  of  the  two  modes  only  at 
threshold  and  above  thresold  with  both  modes  lasing.  In  these  cases 
a change  in  the  photon  density  of  one  mode  will  nave  definite  impact  on 
the  density  of  the  other.  At  threshold,  even  though  the  standard  devia- 
tions are  still  on  the  order  of  the  mean  densities,  this  indicates  the 
onset  of  order  within  the  system.  Below  threshold,  and  above  threshold 
with  only  one  mode  oscillating,  is  very  small  and  the  two  modes  are 

very  nearly  uncorrelated.  Below  threshold,  both  modes  are  dominated  by 
spontaneous  emission.  Therefore,  neglecting  the  small  amount  of  stim- 
ulated emission  that  may  occur,  the  photons  in  one  node  will  have  no 
effect  on  the  photon  number  in  the  other  mode.  Likewise,  above  threshold 
with  one  inhibiting  the  other,  the  inhibited  mode  is  inhabited  almost 
entirely  by  spontaneous  photons  and  can  do  very  little  to  affect  the 
number  density  of  the  lasing  mode. 

Dynamic  Analysis 

General.  Representative  plots  showing  the  evolution  of  the  mean 
photon  number  densities  and  of  the  standard  deviations  and  correlation 
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coefficient  are  shown  in  Figs.  6-13,  Each  figure  comprises  four  plots  i 


nj  vs. 

There  were  no  significant  differences  in  the  behavior  of  the  system  be- 
tween the  two  levels  of  above  threshold  pump  considered  in  this  study. 
Therefore,  only  one  of  these  pump  levels  is  considered  in  this  section. 

The  same  was  true  for  the  two  levels  of  pump  below  threshold.  In  addi- 
tion, the  degree  of  coupling  had  no  real  noticeable  effect  below  thres- 
hold on  the  evolution  of  the  variables  and  caused  only  minor  differences 
at  threshold.  Therefore,  only  one  threshold  case  and 

two  below  threshold  cases  (corresponding  to  matched  and  unmatched  pumping) 
are  considered.  The  remainder  of  the  plots  are  included  in  Appendix  C. 

Below  Threshold.  Figs.  6 and  7 illustrate  the  evolution  of  the  first 
and  second  moments  for  the  two  below  threshold  cases  A41  and  A42,  respect- 
ively. As  predicted,  the  mean  photon  densities  follow  deterministic 
paths.  In  both  cases  the  modes  are  dominated  by  spontaneous  emission. 
Therefore,  they  exhibit  characteristics  of  white  (incoherent)  light. 

The  standard  deviations  s^  and  s^  are  on  the  order  of  the  mean  densities 
throughout  the  transitions  and,  therefore,  experience  the  same  general 
behavior  as  the  mean  photon  densities.  Also,  the  correlation  coefficient 
is  very  near  zero  throughout  the  regions  of  interest  indicating  that  the 
modes  have  little  on  each  other. 

At  Threshold.  At  threshold,  only  the  weakly-coupled  case  A31  has 
real-valued  saddle  points.  The  strongly-coupled  system  A33  has  two 
complex-valued  singularities,  while  the  neutrally-coupled  system  A35 
has  only  two  nodes,  one  stable  and  one  unstable.  Comparing  these  three 
cases,  it  is  seen  that  there  is  little  impact  caused  by  the  presence 
(or  absence)  of  saddle  points. 


n2  t n.  vs.  n2  vs.  s1  j n^  vs.  n2  vs.  s2  ; and  n^^  vs.  ''S.^2 
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Considering  the  case  A35  (Fig.  8),  the  same  general  relationships 
exist  between  the  mean  photon  densities  and  the  standard  deviations  as 
exist  below  threshold.  The  modes  are  still  dominated  primarily  by  spon- 
taneous emission  and,  therefore,  and  s^  remain  on  the  order  of  n^ 
and  . However,  the  correlation  coefficient  no  longer  is  negligible. 
When  either  mode  is  operating  well  below  equilibrium  the  correlation 
between  modes  is  very  small.  As  the  mode(s)  grow(s),  the  effects  of 
saturation  and  coupling  cause  the  modes  to  affct  each  other  and  thus 
they  become  correlated.  Because  of  the  simulated  Q-switch  used  to  deter- 
mine the  initial  variances  and  covariances,  the  curves  associated  with 
large  initial  photon  densities  in  both  modes  start  out  indicating  signif- 
icant correlation  between  modes  which  then  decrease  or  increase  monoton- 
ically  towards  the  equilibrium  value. 

Above  Threshold.  Above  threshold,  the  behavior  of  the  system  is 
strongly  affected  by  the  singular  points  (Table  II).  Consider  first 
the  weakly-coupled  case  A21  as  shown  in  Fig.  9.  Comparing  Fig.  9a, 
the  plot  of  n^  vs.  n^  , to  Fig  1,  its  equivalent  plot  resulting  from 
Lamb's  semiclass leal  treatment,  favorable  agreement  is  noted  with  an 
important  addition.  Fig,  9a  includes  phase  curves  from  points  with  one 
or  both  nodes  in  a photon  vacuum.  The  curves  between  the  saddle  points 

Q Q 

at  approximately  (0,4x10  ) and  (4x10  ,0)  and  the  stable  equilibrium  act 
as  "limit  lines"  of  the  system  to  which  other  phase  curves  are  pulled 
by  the  saddle  points.  This  is  especially  evident  for  those  lines  be- 
ginning with  one  of  the  modes  in  a photon  vacuum  (such  as  those  with  ini- 
tial points  Identified  on  the  n2  axis  by  several  different  characters). 
These  curves  move  quickly  towards  the  closest  saddle  point  before  the 
non-lasing  mode  can  build  significantly,  as  there  are  not  enough  excited 


atoms  to  maintain  the  system  at  Its  current  level.  They  then  follow  very 
closely  the  "limit  line"  in  the  transition  to  stable  equilibrium.  Figs. 

9b  and  9c  Indicate  that  the  occurence  of  anomalous  fluctuations  is 
associated  with  buildup  in  one  or  both  modes  from  a photon  vacuum,  and 
that  the  maximum  fluctuation  occurs  approximately  half-way  along  the  curve 
between  the  stable  and  unstable  equilibria.  In  buildup  from  (n^  , n^) 
-(0,0)  both  modes  experience  the  same  degree  of  fluctuation  as  both 
modes  are  expected  to  grow  at  the  same  rate.  When  only  one  mode  is 
building  from  zero  both  modes  again  experience  the  anomalously  large 
fluctuations,  even  though  the  other  mode  is  already  oscillating.  How- 
ever, the  growing  mode  experiences  a larger  fluctuation.  Comparing 
the  curves  for  s^  and  s,  to  their  associated  curves  for  p in  Fig.  9<i, 
it  is  seen  that  movement  along  the  "limit  line"  is  associated  with  a high 
correlation  between  modes  ( p^ Also,  the  phase  curve  assoc- 
iated with  buildup  from  a photon  vacuum  in  both  modes  moves  towards 
before  relaxing  to  the  equilibrium  value  of  -0.5. 

Case  A22,  representing  weakly-coupled  modes  with  mode  2 inhibited 
at  equilibrium,  is  shown  in  Fig.  10.  As  with  case  A21,  a comparison  of 
Fig.  10a  with  its  equivalent  from  Lamb's  theory  (Fig.  2)  indicates  fav- 
orable agreement  in  n^  and  n^  . There  is  now  only  one  saddle  point,  at 

Q 

(n^  , n2) «(0  , 10  ),  which  affects  the  system  and  all  phase  curves  move 
toward  the  one  equilibrium  where  only  mode  1 lases.  The  anomalously 
large  fluctuations  are  again  associated  with  buildup  from  a photon 
vacuum.  The  slngle-to-dual-mode  transition  of  Fig.  9 is  replaced  by 
transition  from  unstable  to  stable  mode  oscillation,  and  again  both 
modes  experience  the  anomalous  fluctuations,  as  seen  in  Figs.  10b  and 
10c.  However,  only  mode  1 experiences  anomalous  fluctuations  in  build - 
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up  fron  (0,0)  since  it  is  the  only  mode  which  oscillates  at  equilibrium. 
Also,  movement  along  or  close  to  the  "limit  line"  is  again  associated 
with  strong  correlation  between  the  modes,  as  illustrated  in  Fig.  lOd, 

The  case  of  strongly-coupled  modes  pumped  equally,  as  shown  in 
Fig.  11,  is  unique  from  the  other  cases  studied  because  it  has  two 
different  equilibria  which  are  stable  nodes  and  to  which  the  system 
can  relax.  It  also  has  an  unstable  equilibrium  lying  in  the  physical 
region  of  the  (n^,n2)  plane.  This  internal  saddle  point  determines  the 
stable  point  to  which  the  system  will  move.  As  before,  the  evolution  of 
the  mean  photon  densities,  Fig.  11a,  compares  favorably  with  Lamb's 
results  as  plotted  in  Fig.  3.  Since  both  modes  can  be  stable,  depend- 
ing on  one's  current  position,  the  system  does  not  have  an  unstable 
mode  from  which  to  transition,  as  in  case  a22.  Similarly,  it  cannot 
move  to  two-mode  oscillation.  Therefore,  one  might  expect  that  for  this 
case  the  only  anomalous  fluctuations  are  associated  with  buildup  in  one 
of  the  modes  from  a photon  vacuum  while  the  other  mode  remains  essential- 
ly empty.  However,  while  photon  buildup  transitions  do  experience  the 
anomalous  fluctuations,  the  strongest  influence  is  the  internal  saddle 
point.  As  is  evident  in  Figs,  lib  and  11c,  the  phase  curves  experiencing 
the  largest  fluctuations  are  those  passing  closest  to  this  saddle  point. 
Again,  it  is  interesting  to  note  that  the  maximum  fluctuations  occur  not 
nearest  the  unstable  equilibrium,  but  approximately  half-way  between  the 
unstable  point  and  the  stable  point.  As  in  all  the  above-threshold 
cases  studied,  movement  along  the  "limit  lines"  extending  between  the 
stable  and  unstable  equilibria  was  associated  with  high  correlation 
between  the  modes,  or  ^2*-1.0  . However,  for  strongly-coupled  modes, 
due  to  the  stronger  mode  interaction,  this  same  correlation  exists 
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whenever  both  modes  are  oscillating,  as  seen  in  Fig.  lid. 


When  the  strongly-coupled  modes  are  unequally  pumped  such  that  only 
one  mode  is  inhibited,  as  in  case  A24,  the  evolution  of  the  moments  is 
as  shown  in  Fig.  12.  There  is  again  only  one  stable  point  to  which 
the  system  can  relax  and  there  is  only  one  saddle  point  affecting  the 
system.  As  in  case  A22,  anomalous  fluctuations  occur  during  both  trans- 
ition from  unstable  to  stable  mode  oscillation  and  buildup  of  mode  1 
from  a photon  vacuum  (0,0).  However,  the  largest  fluctuations  are 
associated  with  phase  curves  which  begin  with  both  modes  oscillating 
and  then  pass  near  the  saddle  point.  Both  modes  oscillating  under  strong 
coupling  is  a highly  unstable  configuration  and,  therefore,  the  curves 
associated  with  these  points  begin  with  large  uncertainties  on  the  order 
of  the  mean  photon  densities.  As  the  system  evolves  toward  the  stable 
equilibrium,  S£  maintains  the  same  relative  value  throughout  the  trans- 
ition (Fig.  12c).  Likewise,  for  those  curves  which  approach  the  "limit 
line"  before  the  half-way  point  where  the  maximum  fluctuations  are  ex- 
pected to  occur,  s^  stays  on  the  order  of  n^  and  as  n^  grows,  so  does  s^ 
(Fig.  12b),  Once  the  maximum  point  is  reached,  s^  decreases  rapidly  and 
approaches  its  relative  low  value  characteristic  of  stable  lasing  modes. 

Looking  at  Fig.  13,  there  are  only  minor  differences  between  the 
evolution  curves  for  neutrally-coupled  modes  and  those  for  strongly- 
coupled  modes  under  unequal  pumping.  For  the  neutrally-coupled  modes 
(case  A25)  the  effect  of  the  saddle  point  is  not  as  pronounced.  The 
phase  curves  do  not  approach  as  close  to  the  saddle  point  and,  except- 
ing those  which  begin  on  the  axis,  do  not  approach  the  "limit  line" 
as  quickly.  Therefore,  even  though  the  two-mode-oscillation  initial 
points  have  associated  uncertainties  on  the  order  of  the  photon  densities, 
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the  loops  in  the  curves  of  Fig.  12b  do  not  occur  in  the  s^  curves  of 
Fig.  13b. 

It  is  interesting  to  look  at  the  general  shape  of  the  "limit  lines" 
relative  to  the  mode  coupling.  When  the  modes  are  weakly  coupled  and 
and  the  system  is  in  the  vicinity  of  a saddle  point,  one  of  the  mode  den- 
sities will  be  near  zero.  As  this  mode  grows,  the  other  mode  density  will 
decrease  as  both  modes  pull  towards  their  own  respective  equilibria. 
However,  the  decrease  required  will  be  less  than  the  growth  of  the  first 
mode.  Therefore,  the  "limit  lines"  are  concave  downward  or  towards  the 
origin.  On  the  other  hand,  the  "limit  lines"  are  concave  upward  for 
strongly-coupled  modes  as  the  opposite  relative  responses  between  modes 
are  encountered.  'When  the  system  is  neutrally-coupled,  the  "limit  line" 
becomes  a straight  line  showing  no  preferred  concavity,  just  as  the  system 
shows  no  preferred  method  of  saturation. 

In  comparing  Figs.  9a,  10a,  and  11a  to  the  results  of  Lamb's  theory 
(Figs,  1,  2,  and  3)  the  one  important  difference  noted  was  the  inclusion 
of  phase  curves  originating  from  a photon  vacuum.  There  is  one  other 
significant  difference  between  these  two  models.  Lamb's  model  pre- 
dicted that  any  point  along  lines  1^  and  12  (Eqs  (43)  and  (44))  would  be 
stable  with  deviations  away  from  either  line  decaying  back  towards  that 
line.  However,  when  the  spontaneous  emission  is  included,  this  stability 
is  destroyed  and,  as  seen  in  Fig.  13,  the  system  moves  to  equilibrium 
at  a single  point. 

Mode  Buildup  from  Vacuum.  Of  particular  interest  is  the  buildup 
of  the  mode  densities  from  a photon  vacuum.  Figs,  14a,  14b,  and  14c 
illustrate  this  buildup  in  time  for  several  levels  of  pump  and  degrees 
of  coupling.  Fig.  14c  shows  the  buildup  in  both  modes  from  a photon 
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vacuum  for  the  five  different  above  threshold  cases  shown  in  Figs.  9-13. 
Fig.  14b  corresponds  to  Fig  8,  and  Fig.  14a  does  the  same  for  the  cases 


of  Figs.  6 and  7,  Below  threshold,  the  modes  grow  independently  with 
the  time  relationship 

n;  = «in  ( 1 - e"** ) (74) 

where  n^  is  the  equilibrium  density  for  that  mode  and  a^  is  a decay 
constant.  This  decay  constant  is  equal  to  the  net  losses  of  the  mode, 
or  the  difference  between  the  loss  parameter  h<  and  the  stimulated  gain 
parameter  g£^  . Coupling  and  saturation  have  negligible  effects.  At 
threshold,  according  to  Eq  (74)  both  modes  should  increase  linearly  with 
no  damping  or  aplification  effects,  since  the  stimulated  gains  equal 
the  losses.  This  does  occur  initially.  However,  as  the  mode  densities 
continue  to  increase,  saturation  sets  in  and  the  rate  of  buildup  de- 
clines as  equilibrium  is  reached.  Above  threshold,  both  modes  exper- 
ience initial  exponential  rises  with  growth  constants  equal  to  the  net 
gain  of  the  system.  As  they  approach  and  pass  the  threshold  equilibria, 
saturation  and  coupling  begin  to  modify  this  rise,  gradually  slowing  the 
rate  of  increase  and  even  causing  a decrease  for  the  mode  which  is  in- 
hibited from  lasing  at  equilibrium. 

The  time  to  equilibrium  is  primarily  determined  by  the  size  of  the 
net  gain  or  loss  of  the  system.  The  closer  the  system  is  to  threshold 
operation,  the  longer  it  will  take  for  the  system  to  buildup  to  equil- 
ibrium. Below  threshold,  the  net  loss  (or  the  decay  constant)  js  the  only 
determining  factor,  whereas  above  threshold,  and  even  at  threshold,  the 
effects  of  coupling  and  inhibition  speed  up  or  slow  down  the  process. 
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For  example,  consider  the  times  to  equilibria  as  shown  in  Figs.  14a- 
14c.  At  threshold  (Fig.  14b) , this  time  was  on  the  order  of  1 msec. 

Below  threshold  (Fig  14a),  it  took  about  jj^sec  for  the  system  to  reach 
an  equilibrium.  In  Fig.  14c  the  times  to  equilibrium  were  in  the  range 
of  5-lOyu.sec, 

As  with  the  photon  densities,  the  buildup  of  the  mode  densities 
is  used  to  illustrate  the  evolution  of  the  standard  deviations  in  time. 

See  Fig.  14d.  Since  the  standard  deviations  for  systems  at  and  below 
threshold  are  very  close  to  the  mean  photon  densities,  only  the  above 
threshold  cases  represented  in  Fig  14c  are  considered.  Also,  the  strongly- 
coupi>  equally-pumped  system  experienced  fluctuations  two  orders  of 
magnitude  greater  than  the  fluctuations  for  the  other  four  cases  in- 
volved. Therefore,  the  time  behavior  of  its  standard  deviations  is  not 
included  in  Fig,  14d,  Comparing  the  curves  of  Fig.  14d  to  Fig  14c,  it  is 
seen  that  the  fluctuations  do  peak  approx ^nately  half-way  to  equilibrium, 
as  predicted  and  seen  before,  with  the  case  of  two-mode  oscillation 
peaking  farther  towards  the  equilibrium  than  the  cases  of  single-mode 
oscillation.  The  height  of  the  peak  Is  dependent  primarily  on  the  pump- 
ing, as  different  pump  parameters  produce  different  photon  densities. 

There  is  also  a dependence  on  the  mode  coupling  for  the  same  pump  level. 

For  each  of  the  curves  in  Fig  I4d,  the  peak  represents  approximately 
the  same  relative  standard  deviation  ( s^n^  ).  Therefore,  since  the 
peak  occurs  when  the  photon  density  is  about  half  that  of  the  equilibrium 
value,  the  heights  of  the  peaks  will  have  the  same  relative  values  as  do 
the  equilibrium  photon  densities. 
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V Conclusions  and  Recommendations 


Conclusions 

Based  on  the  results  of  this  study,  several  conclusions  can  be 
drawn.  The  most  important  of  these,  without  which  the  rest  would  be 
meaningless,  is  that  the  method  used  in  this  study  of  photon  statistics 
is  valid.  The  lasing  process  can  be  described  adequately  by  the 
Chapman-Kolmogorov  master  equation  and  the  proposed  transition  rates 
used  in  solving  the  master  equation  are  appropriate.  Accepting  this, 
the  other  conclusions  are  i 

1) .  Any  two-mode  laser  system  has  at  most  two  stable  equilibria, 
with  all  but  some  strongly -coupled  systems  having  only  one.  This  in- 
cludes systems  with  neutrally-coupled  modes,  contrary  to  the  predictions 
of  the  Lamb  model. 

2) .  The  laser  system  experiences  large,  anomalous  fluctuations  in 
one  or  both  modes  in  three  specific  types  of  transition.  First,  when- 
ever the  system  builds  up  from  a near  photon  vacuum  such  as  at  initial 
startup,  the  mode  or  modes  which  will  be  oscillating  at  equilibrium  will 
experience  the  anomalous  fluctuations.  Second,  whenever  only  one  mode  is 
oscillating  and  both  modes  oscillate  at  equilibrium,  both  modes  will 
experience  the  anomalous  fluctuations  during  the  transition  to  equil- 
ibrium. Third,  whenever  an  unstable  mode  is  oscillating  alone  both  modes 
will  have  anomalous  fluctuations  in  the  transition  from  the  unstable 

to  the  8 table  mode. 

3) .  Strongly-coupled  laser  systems  operating  with  two  modes  os- 
cillating experience  large  fluctuations,  especially  when  there  are  two 
possible  stable  points  to  which  the  laser  system  can  transition. 
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4) .  A mode  experiences  normal  fluctuations  when  the  system  is 
below  threshold,  when  the  mode  both  begins  and  ends  the  transition  below 
threshold,  and  when  one  mode  ceases  lasing  while  the  other  mode  con- 
tinues to  lase  (excepting  the  case  where  two  physical  stable  points 
exist) . 

5)  . Inhibited  modes  above  threshold  experience  the  same  type  of 
fluctuation  as  if  they  were  below  threshold. 

6) .  The  anomalous  fluctuations  are  accompanied  by  a strong  linear 

relationship  between  the  mode  densities  ( ~ 1 )• 

7) .  There  is  significant  correlation  at  equilibrium  between  the 
nodes  only  at  threshold  and  above  threshold  with  both  modes  oscillating. 

Recommendations  for  Further  Study 

1) .  A better  qualitative  understanding  of  the  effect  of  coupling 
could  be  gained  by  solving  several  different  cases  for  various  values  of 
the  coupling  parameters. 

2) .  Operation  very  close  to  threshold  should  be  given  a closer 
look,  possibly  resulting  in  a better  definition  of  "near  threshold". 

3) .  This  study  was  limited  to  constant  pump  parameters  which 
turned  on  or  changed  instantaneously.  Therefore,  a similar  analysis 
could  be  done  allowing  the  pump  parameters  to  be  functions  of  time. 

This  could  then  allow  for  analysis  of  pulsed  operation  as  well  as 
continuous  wave. 

^0  • This  analysis  could  be  repeated  for  other  proposed  expressions 
for  the  transition  rates, 

5).  The  equation  development  used  in  this  study  should  be  ex- 
panded to  include  multimode  laser  operation  and  the  general  results 
for  two-mode  operation  extended  as  appropriate. 


59 


Bibliography 


1.  Arecchi,  F.  T.  and  A.  M.  Ricca.  "Statistical  Properties  of  a Many- 
Mode  Laser,"  Physical  Review  A,  l£i  308-318  (January  1977). 

2.  Arecchi,  F.  T.  et  al.  "Time-Dependent  Statistical  Properties  of 
the  Laser  Radiation,"  Physical  Review  Letters,  l£i  1168-1171 
(13  November  1967). 

3.  Haken,  H,  "Cooperative  Phenomena  in  Systems  Far  From  Thermal  Equil- 
ibrium and  in  Nonphysical  Systems,"  Reviews  of  Modern  Physics,  47 > 
67-121  (January  1975). 

4.  Kubo,  R.  "Non-linear  Relaxation,"  Statistical  Physics—  S tatPhys 
13.  Annals  of  the  Israel  Physical  Society,  Vol  2.  Part  1_,  edited 
by  D.  Cahib,  et  al.  New  Yorki  Adam  Hilger,  Bristol  and  the  Israel 
Physical  Society  in  association  with  the  American  Physical  Society, 
1978. 

5.  — - — . "Relaxation  and  Fluctuation  of  Macrovariables,"  Synergetics  - 
Cooperative  Phenomena  in  Multi-component  Systems,  edited  by  H.  Haken. 
Stuttgarti  Teubner,  1973. 

6.  Kubo,  R.  et  al.  "Fluctuation  and  Relaxation  of  Macrovariables," 
Journal  of  Statistical  Physics.  2*  51 "96  (September  1973). 

7.  Lamb,  W,  E.,  Jr.  "Theory  of  an  Optical  Maser,"  Physical  Review, 

134;  A1429-A1450  (15  June  1964). 

8.  Mandel,  P,  "Quantum  Multi-mode  Laser  Theory,"  Physlca  B & C. 
Netherlands.  82»  353-367  (April  1976) . 

9.  Risken,  H.  "Statistical  Properties  of  Laser  Light,"  Progress  in 
Optics.  VIII i 241-296  (1970) 

10.  Sargent,  M,  III  et  al,  "Buildup  of  Laser  Oscillations  from  Quantum 
Noise,"  Applied  Optics.  2423-242?  (November  1970). 

11.  — - — . Laser  Physics.  Reading,  Massachusetts!  Addison-Wes-ley 

Publishing  Co.,  197^.  * 

12.  S legman,  A.  E.  An  Introduction  to  Lasers  and  Masers.  New  York* 
McGraw-Hill  Book  Co.,  1971. 

13.  Skappel,  J.t  Boeing  Aircraft.  "FORTRAN  IV  Subroutine  BLCKDQ," 
Coordination  Sheet  No.  PS-769,  22  March  1968. 

14.  Standard  Math  Tables  (17th  Edition),  edited  by  Samuel  M.  Selby. 
Cleveland,  Ohioi  The  Chemical  Rubber  Co.,  1969. 


60 


15.  Stoker,  J.  J.  Nonlinear  Vibrations  In  Mechanical  and  Electrical 
Systems.  New  Yorki  Interscience  Publishers,  Inc.,  1950. 

16.  van  Kanpen,  N,  G.  "A  Power  Series  Expansion  of  the  Master 
Equation,"  Canadian  Journal  of  Physics.  39 » 551-567  (1961). 


1 


61 


l 


Appendix  A Derivation  of  the  Evolution  Equations 


Let  X ■ be  a vector  variable  representing  the  photon 

numbers  of  a two-mode  laser  system.  It  is  assumed  that  X is  a stochas- 
tic process  as  it  evolves  in  time.  Further,  X is  assumed  to  be  Marko- 
vian and,  therefore,  can  be  described  by  the  Chapman-Kolmogorov  equation 
in  the  form 

= -/w(x,r>(X,t)<lr  + /w(X-rir>('X-C,i)<k  (7?) 

where  P(X,t)  is  the  probability  distribution  of  X at  time  t and  W(X,r) 
is  the  transition  rate  for  a jump  from  X by  an  amount  r.  Since  the 
values  of  X will  be  discrete,  the  integrations  are  understood  as 
summations . 

X,  as  defined  above,  represents  the  extensive  nacrovariables  of 
the  laser  system,  each  of  which  is  made  up  of  contributions  from  many 
elementary  units  - in  this  case,  photons.  However,  it  is  assumed  that 
the  transition  rate  W(X,r)  is  determined  by  the  internal  state  of  the 
system,  rather  than  by  X.  Therefore,  the  transition  rate  takes  the 
form 

(7 C,) 

where  _f L is  the  size  of  the  system  and  x is  the  normalized  intensive 
macrovariable  vector  corresponding  to  X and  defined  by 

a = X//l  (ti) 
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With  this  assumption,  the  master  equation  (75) • rewritten  for 

, becomes 

JTvt  4-Jw(i-X.^r')p(x-c_|t)dr  (ts) 

By  expanding  the  second  integral  on  the  right-hand  side  of  Eq  (78)  in  a 
power  series  about  x,  Eq  (78)  can  be  written  formally  as 

iiWfti,*)*  ~f  I w(s-.c>Ki‘,*')J.c  (rO 

or,  maintaining  the  expanded  form 

(ir)  Tij/ <k(£  • ?)"<-(<, c)rt&,t)  o<o 

where  £ is  the  transpose  of  r. 

Following  the  work  of  Kubo  (Ref  5«29)»  a solution  of  Eq  (80)  of 
the  form 

p(xit')  = C ocpCjlft's.t)]  O') 

is  assumed.  Then,  the  time  derivative  of  P(x»  t)  is 

= \t  [c  l/^M]  = ji  ,b)  ( k) 

and  the  integrand  on  the  right-hand  side  of  Eq  (80)  becomes 
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(c_-  v)  = (f.*  v)  uj(y  C -a-^pC/2--f(x)tS)] 

=Jln *vf) 

+■  0(-An"')  ^ 

o 

Therefore,  to  order  _TL  , Sq  (80)  becomes 

= £ ('0*4r / ■?()"  (h) 

n-i  • ' 

and  (79)  becomes 

TJ  ((*,*)  = -/ dc^.c)  [l  - W-r  • 2-f(s,t))]  (Kf) 


To  solve  this  equation,  a solution  of  the  form 


-f(*,t)  = 8(2)*)  , 


00 


is  assumed,  where  Jf(t)  is  the  expected  value  of  x and  ;s  is  the  deviation 
of  x away  from  This,  then,  transforms  (85)  into 

= ~ dr  ^(^-*-2  >r)  |j  -^pfe-v-j)]  ($7) 

Assuming  4 (*T)|  • the  transition  rate  w(yfa,r)  can  be  written  as 


(ts) 


which  then  transforms  (87)  into 


~ J ^S.  "Z*  Yju^(^^rV)  [j  - £.*  V'j'j]  ( 


Expanding  the  exponential  terras,  Eq  (89)  becomes 

+ £ • v,  J eYl  • r 


(lo) 


Noting  that  the  first  moment  of  the  transition,^*d£oj(^r)_r  f is  just 
the  average  velocity  of  £ when  the  system  is  in  the  state  £,  and  that  jr(t) 
is  already  an  average  or  expected  value,  it  is  easily  seen  that  j(t) 
satisfies  the  condition 

^ (V)  - J d_r  (^0 

Therefore,  the  terras  in  Eq  (90)  linear  in  Vj^vi^O  vanish  and  the 
resultant  equation  is 

7T  s(i,t)  (dr  wt^,c)(r  •vlS)'1 

+ (<u) 

n»i  ^ ^ 
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Assuming  an  approximate  gausslan  solution  of 

<3  m = _ x i ■ a (o  • i + od  ir) 

where  A Is  the  Inverse  of  the  covariance  tensor,  then 

£.  AW  * £)  + 0(\  ^1  ) 

= \i^i)  +*  0(li.O 

- vA3*  "** 

However 

A = A 0*) 

/ ' ^ 1 y * 

and,  therefore 

r-vJ'j(?,t)=  - Ci.A.^5^  +■  0(lil  ) 

= -f-A-a  + 0(U|3)  00 

Inserting  (93)  and  (96)  Into  Eq  (92)  and  taking  the  leading  terms  to 
order  (i |*  , Eq  (92)  becomes 

-i_£.  A*  2,  = y|*A  ' Jdr  fj’A'j. 

+ 1 • A * ^ (0|7^ 
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This  leads  directly  to 


A/t)=  -A  (is) 

which  can  then  be  made  symmetric  (recovering  the  symmetry  lost  by  using 
(95)  to  combine  the  terms  of  Eq  (9^)  by  noting  again  the  equivalence 
of  Eq  (95).  Therefore 

VjJdrw(i,r')r -A  >^)  ^ 

and  Eq  (98)  becomes 

Aw  * "A  'J dr  «*»(},  L^cc-A  -vjdr  uj(^Or  . /\ 

~ A • ( Wn)  (100) 

An  equivalent  expression  for  the  covariance  tensor  B(t)  is  obtained 
by  first  noting  several  tensor  identities.  By  definition 

S “ A_  (/oO 

and,  therefore 


B = -BAB 


(|o4) 


Defining  c ^ (^)  and  c 2 (jr)  as  the  first  and  second  moments  of  the 
transition  rate  w(^fr) 


c,m 


ol) ( y t r)  r 

(tofr) 

w(<4  ££ 

(l  0 O 

and  inserting  (100)  into  Eq  (104),  the  equation  for  the  covariance 
tensor  is 


bW  = tx(^)  + 


(10^ 


or 

yrA./5)  +^c'“-  W BV(S  (10 o 


To  obtain  the  covariance  tensor  for  x,  (107)  must  be  divided  by 
This  is  seen  by  first  inserting  (93)  into  (87)  to  get  (keeping  only 
the  leading  term  in  | i ) z ) 


P ( * } ~ C 2.  2.  • A • 2] 

- C jucp  []” ■£  i • (-/I^V  (toO 
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Letting 


cr~'(0  - ~fl  A (no) 

be  the  inverse  covariance  tensor  for  x,  the  covariance  tensor  becomes 

i(t)  = (in) 

Similarly 

X(t)  = ji  (m) 

where  AW  is  the  covariance  tensor  for  X.  Therefore,  by  substituting 
(111)  back  into  Eq  (107)  the  equation  for  the  covariance  tensor  for 
x results  < 


iO)  + +(2,£,(a.))sr  (ns) 

where  c^  and  c2are  as  defined  in  Eqs  (105)  and  (106). 

Thus,  Eq  (113)  and  Eq  (9?)»  which  can  be  rewritten  as 

£W  = c.,(s)  (|I4) 

are  the  general  equations  governing  the  evolution  of  the  first  and 
second  moments  of  the  normalized  intensive  macrovariable  x.  Equation 
(114)  defines  the  most  probable  path  y(t)  which  the  variable  x will 
follow,  while  Eq  (113)  determines  a measure  of  the  fluctuatias  of  the 
variable  around  this  average  path. 


Appendix  B Linearization  of  the  Photon  Rate  Equations 


To  linearize  the  photon  rate  equations  (36)  and  (37) t consider  a 
snail  disturbance  around  an  equilibrium  point.  Let 


and 


a , - n,  +*  ^ 


- < 4- 


(l/O 


CmO 


where  n°  and  n|  are  the  equilibrium  point  values.  Substituting  these 
into  Eq  (36),  the  rate  equation  becomes 


* *°  , ’ ' - Si.  4 (ft*  ) - \\  ( 

n''-n'+n'_  ' 


:)  (in) 


5..  4 A j -4-  A, 


1 +•  4<x<\ 


(l\i) 


Now 


x + A X 

Therefore,  letting 


* ~ A/ 

X1 


X = l 4 ©<ufl*  4-  Xo.Oi0 


and 


AX-  <*.,  A,  4-  ©<,2. 
Eq  (118)  becomes 


(in) 


Ozo) 

(.1  2-0 


• 1 1 

( + <Lu’A.°'*“SuAt  . 

l(l  4^0°  - «><«  _ 1 

h *i°-  a 

p 

>1 

1 

( 1 4 A * 4 °^vj.  A j_  ^ 

•i'u  "|m» 
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Separating  out  zero  order  terns 

(^.>+  -<A,' -o(uK 

( * 4 + <*,t  <a“^ 


-vt^r 


\[  (i4°w°v<*>t.ft  __  ViX 

( 1 +■  <x  ox  +- 


(1^3) 


= f _h,n; 

\ 1 4c4,vA,  4^u.Al 


! 4 ^iifli  4o<'ia.0t. 

| -f-  4 <XtL  'A* 


, q^(l^v,A^,,.nl-*„rt,'-**AO  _ ^,a?  4°4iQz)  - K ,0/  (424) 

( i +**:  i+^n.° 

Now,  the  first  bracketed  factor  of  the  first  term  in  Eq  (124)  is  just 
which  is  equal  to  zero.  Therefore,  the  entire  first  term  vanishes  and 


the  rate  equation  becomes 


( ( 4<n“  4 °( u.nt  ^ 


- w (us-) 

I -J-tfufl,0  +■  C^,t  A° 


Since  the  disturbance  was  originally  assumed  to  be  small,  only  the  first 
order  terms  are  of  significance;  therefore,  separating  out  all  higher 


order  terms,  Eq  (125)  becomes 


H.O.T.  XO 


* ' _ Sufi,  t-oCn.A*\  __  3z\Av  1 

' ( l 4 4 ( I 4 °Ai.  A*.  ^ 

V^A,  (°A>AV  ^ ^nAi.)  i_.  Vl  j A ^ 

1 4 <*,,0,° 


GzO 


71 


» 


Appendix  C Additional  Evolution  Curves 
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